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Abstract 



The paper derives an equation for the Cauchy transform of the solution of a 
free stochastic differential equation (SDE). This new equation is used to solve 
several particular examples of free SDEs. 



1. Introduction 

Free stochastic differential equations generalize classical stochastic differential equations 
to the setting of free probability. Here is an example of such an equation: 

X t = a{X t )dt + b*(X t )(dZ t )b(X t ) 

In this equation, Xt is a self-adjoint operator, a(Xt) and b(Xt) are operator-valued functions 
of Xt, and the driving noise Z t is an operator process with free increments. That is, the incre- 
ments Z s — Z t , s > t, are assumed to be free from past realizations of Z t . The process Z t is 
usually the free Brownian process, in which case the increments have semicircle distributions; 
however, other choices are possible. 

Informally, the reader may think about X t as very large random matrices and Z t as matri- 
ces with independent Gaussian random variables as entries. These entries follow independent 
Brownian motions and we are interested in the law of the eigenvalues of X t . The free proba- 
bility theory is a convenient abstraction which intends to model the situation when the size of 
the matrices is very large. 

The study of free stochastic differential equations ("free stochastic calculus") is more dif- 
ficult than in the classical case because of non-commutativity of coefficients and noise. This 
paper contributes by developing a new tool for the analysis of these equations. 

The idea of free stochastic calculus was first suggested in (Sj. It was later developed and 
formalized in [KS], [BJ, [BSa], and [AJ, which introduced stochastic integration with respect 
to free Brownian motion as a rigorous basis for free stochastic calculus. They also derived an 
analog of the Ito formula, which allows us to obtain identities like the following: 

( a [W?(dW t ) + W t {dW t )Wt + (dW t )W?] =W z a -2 [ W t dt, 
Jo Jo 
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where Wt denotes the free Brownian motion (the Wigner process). An analogous formula in 
the classical situation is 



/ 3B 2 (dB t ) = B 3 a -3 B t dt, 
Jo Jo 



where B t is the standard Brownian motion. Note the different coefficient before the integral 
on the right-hand side. 

The classical Ito formula is very helpful in the study of stochastic differential equations. 
Unfortunately, the range of applicability of the free Ito formula is smaller. This difficulty 
calls for a different method applicable to those free SDEs, which are not solvable with the Ito 
formula. One possibility is to seek an equation for the evolution of the spectral distribution of 
the solution. Such an equation was derived by Biane and Speicher in MBSbl for the equation 

X t = a(X t )dt + dW t . (1) 

They showed that the density of the spectral probability measure of X t , which we denote p t 
and which generalizes the eigenvalue distribution of a matrix, satisfies the free Fokker-Planck 
equation: 

Here H denotes a multiple of the Hilbert transform: 

Hu(x) := p. v. / — ^—dy. (3) 
J x-y 

Still, the approach through the free Fokker-Planck equation has its own disadvantages. 
First, it is applicable only to equations that have the special form (QQ), that is, only to equations 
with the constant diffusion coefficient. Second, the free Fokker-Planck equation (0) is not a 
bona fide partial differential equation since it includes the Hilbert transform operator. For this 
reason, it is somewhat difficult to solve this equation. 

The purpose of this paper is to approach the free stochastic equations by deriving a differ- 
ential equation for the Cauchy transform of the solution. 

Recall that the resolvent of operator X t is defined as the operator-valued function of a 
complex parameter Gt (z) := (Xt — . The Cauchy transform of X t is defined as the 
expectation of the resolvent: gt(z) := E[(X t — z)~ ]. It is useful because the knowledge 
of the Cauchy transform is sufficient to recover all properties of the spectral probability 
distribution of X t . It turns out that if X t solves 

dX t = a(X t )dt + b(X t )(dW t )c(X t ), 

then gt(z) satisfies the following equation: 

^ = -E(a t G 2 t ) + E(b t c t G t )E(b t c t G 2 t ), (4) 

where we use at, bt, and q to denote a(X t ), b(X t ), and c(X t ), respectively. This is the 
statement of Theorem |3.2| below. 

Equation (@]) is not a usual differential equation since it involves expectations. In general, 
these expectations are difficult to compute because the coefficients at, b t , and ct, and the 
resolvent Gt are not free from each other. However, if the coefficients are polynomials, it is 
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possible to perform further reduction to a differential equation as we will show in Proposition 

m 

As we just said, the knowledge of the Cauchy transform can be used to recover the spectral 
probability distribution. In particular, the free Fokker-Planck equation © can be derived from 
dU) as will be shown in Corollary 13.61 

In certain cases it is not possible to compute the Cauchy transform explicitly, but it is possi- 
ble to detect the behavior of its singularities. This knowledge can provide us with information 
about the support of the spectral distribution. In particular, it can show us how the norm of the 
solution grows. 

For a simple example of this approach, let us consider the well-known case of the free 
Ornstein-Uhlenbeck equation: 

dX t = -OXtdt + <rdW t . (5) 

For this equation, it is easy to compute the Cauchy transform using equation (@]) and re- 
cover the known result that for positive 9, the spectral probability distribution of X t con- 
verges to a stationary solution, which is a semicircle distribution supported on the interval 

[-*y/W\,vVV\0\]' 

As a more difficult example, consider the equation 

dx t = ex t dt + xt /2 (dw t )xt /2 , 

which can be thought of as a free analog of the equation for the "geometric Brownian motion", 
dxt = Oxtdt + xtdB t . 

Let Xq = I. Equation (01) leads to the following differential equation: 

with the initial condition g(0,z) = (1 — z)^ 1 . The method of characteristics gives us a 
functional equation for the Cauchy transform: 




Figure 1. 6 = 1/2 

While it is difficult to extract an explicit analytical formula for the solution of this equation, 
we can investigate how the support of the distribution changes with time. It turns out that for 
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Figure 4. 9 = — 1, Large times 
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9 < 0, the support of the distribution shrinks to zero. If 9 is between and 1, then the lower 
boundary of the support decreases to zero and the upper boundary grows exponentially fast to 
infinity. If 9 > 1, then both the lower and upper boundary of the support grow exponentially 
fast to infinity. 

We can solve equation © numerically and recover the density of the spectral distribution 
by using Stieltjes formula. Figures [Q4] show the evolution of the density for various values 
of parameters and illustrate the complexity of the behavior of the spectral distribution. For 
example, Figure |2] shows that even if 9 > 1, the spectral distribution does not approach 
infinity immediately. There is a transition period in which a significant portion of the spectral 
distribution remains below A = 1. Similarly, Figure [3] shows that for 9 < 0, the distribution 
does not collapse to zero immediately. Only when time increases, the distribution begins the 
rapid approach to zero, as shown in Figure [4] 

Let us compare this result with the classical analog. By using the Ito formula, it is easy to 
show that the solution of the classical equation for the geometric Brownian motion is 

x t = exp{(6--)t + B t }. 

Hence, with probability 1, the classical solution will decrease exponentially to zero if 9 < 1/2, 
and will grow exponentially to infinity if 9 > 1/2. However, the support of the solution 
distribution is (0, oo) for all t. This is quite unlike the behavior of the free SDE solution. 

Note that the equation for the geometric Brownian motion can be generalized to the free 
probability setting in a different way: 

dx t = ex t dt + x t dw t + (dw t )x t . 

The behavior of the solution of this equations is quite different. In particular, the ratio of the 
standard deviation to the expectation is y / 2(e 2t — 1). This ratio grows exponentially fast with 
t, quite unlike the previous example, where this ratio equals yi. Unfortunately, the partial 
differential equation associated with equation is more difficult to solve and it is not clear 
whether the solution becomes unbounded in finite time. 
Finally, let us consider the following equation: 

dX t = kX t (dW t )X t , 

and let the initial condition be Xq = al. For this equation it is possible to write an explicit 
formula for the spectral distribution of the solution. An interesting feature of this equation is 
that the solution blows up in finite time r = (ak)~ 2 , by which we means that the operator 
norm of the solution becomes infinite as t approaches r. 

Another interesting feature is that as time t approaches r, the spectral distribution con- 
verges to a fixed distribution. If k = a' 1 , then the density of this distribution is 
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FIGURE 5. Density functions for t = 0.1, 0.2, 0.3, 0.4. 
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Figure 6. Density functions for t = 0.5, 0.6, 0.7, 0.8 

which is supported on the interval [1/4, oo) . Otherwise, it is a scaled version of this distri- 
bution. The behavior of the solution density for various times is illustrated in Figures [5] and 

m 

Several specific classes of free SDE have already been investigated in the literature. Biane 
and Speicher in [BSb] and Gao in (Gl studied the free Ornstein-Uhlenbeck equation ©. Biane 
and Speicher proved that its solution converges to a stationary process with a semicircle dis- 
tribution. Gao considered free Ornstein-Uhlenbeck processes with a free Levy driving noise, 
and showed that every self-decomposable probability measure on the real line can be realized 
as a distribution of such a process. 

Capitaine and Donati-Matin in RCD-M1 defined the free Wishart process and found that it 
satisfies the free SDE of the form: 

dX t = Xdt + y/XtdZt + dZ* t yfx t , 
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where Z t is the complex Wigner process. Demni in jDj studied the so-called free Jacobi 
processes which satisfy equations similar to the following: 

dX t = (91 -X)dt + V/ - XtdZty/Yt + y/xldZfy/l-Xt. 

With exception of the free Ornstein-Uhlenbeck process, we study a different set of free 
SDEs, and we approach these equations with a different point of view based on the differential 
equations for the Cauchy transform. 

For the free Ornstein-Uhlebeck process, our results agree with results in MBSbll . 

The rest of the paper is organized as follows. Section [2] provides preliminary information 
about free stochastic integration and Ito formulas. Section [3] describes main results. . In par- 
ticular, Section [3TT1 is devoted to a local existence and uniqueness result. Section 13^21 presents 
general results about the Cauchy transform of the solution and Section l331 provides examples. 



2. Free Stochastic Integration 

2.1. The free Brownian motion. For the basics of free probability theory we refer to HVDNl 
and [NS]. All operators that we consider belong to a non-commutative W* -probability space 
(A, E), that is, to a von Neumann operator algebra A with a faithful normal trace E. We 
denote the usual operator norm by ||X||, and the L 2 -norm by ||X||2 := \/ E(X*X). 

The spectral probability distribution of a self-adjoint operator X € A is a probability 
measure /ionl such that 

E (x k ) = I x k n (dx) . 
Its Cauchy transform is the function 

9x (z) 



fi (dx) 

R x ~ z 

It can be defined directly in terms of operator X as the expectation of the resolvent: gx (z) = 
E[G X 0)], where G x (z) := (X - z)" 1 . The probability measure fj, can be recovered from 
its Cauchy transform by the Stieltjes inversion formula: 

p, (B) = — lim / Img (x + ie) dx, (7) 

provided that B is Borel and \i (dB) = 0. 

This fact is the starting point of our approach, since we will study the evolution of the 
Cauchy transform as a tool to investigate the evolution of the corresponding probability mea- 
sure. 

The most important concept in free probability theory is that of free independence. Let 
Ai denote an arbitrary element of algebra Ai . The sub-algebras Ai,..., A n of algebra A 
(and operators that generate them) are said to be freely independent ox free, if the following 
condition holds: 

E(A i( x- ) . . . Aj( TO )) = 0, 
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provided that E(A i r s \) = and i(s + 1) ^ i(s) for every s. Two particular consequences of 
this definition is that (i) E(AB) = E(A)E(B) if A and B are free, and (ii) 

E{AX X AX 2 ) = E(A 2 )E(X 1 )E(X 2 ), (8) 

if A is free from X\ and X 2 and E(A) = 0. 

The, free Brownian motion, or the Wigner process, is a family of operators Wt, where t > 0, 
that satisfies the following properties: (1) Wo = 0; (2) the increments of Wt are free in the 
sense of Voiculescu, i.e., if t > s, then Wt — W s is free from the subalgebra W s which is 
generated by all W T with r < s , and (3) the spectral distribution of Wt — W s is semicircle 
with zero expectation and variance t — s. 

The choice of Wt is not unique, and in the rest of the paper we assume that a particular 
realization of Wt is fixed. 

2.2. Free stochastic integral. Ito-style free stochastic integration with respect to the free 
Brownian motion was defined and studied in MKSII and ffBSal . Their results show that under 
certain assumptions on the operator coefficients at and bt , it is possible to define the integral 

/= f at(dW t )b t , 
Jo 

where Wt is the free Brownian motion. 

Let us briefly recall the construction of the integral. For details of the construction, the 
reader is advised to see Definition 2.2.1 and Section 3 in [BSa] and Section 3 and Theorem 14 
in (A]. Suppose that at and b t are functions of W T , t < t. That is, let at and bt belong to the 
sub-algebra Wt. Assume also that max{||at|| , ||&t||} < C for all t £ [0, 1] and that t — > at 
and t — > bt are continuous mappings in the operator norm. Let to, . . . ,t n and t\, ■ ■ ■ ,T n be 
real numbers such that 

o = t <h<...<t n = i, 

and 

< Tfe < tk-1- 

We denote the set of to, ... , t n and n, . . . , r n as A. Let 

d(A) = max (t k -r k ). 

l<k<n 

Consider the sum 

n 

I(A) = Y j a Tl (W u -W u _ 1 )b Ti 
i=l 

It turns out that as d(A) — > 0, the sums 1(A) converge in operator norm and the limit 
does not depend on the choice of ti and 73. The limit is called the free stochastic integral 
and denoted as Jq 1 at(dWt)bt- An important point in the proof of convergence is that the 
convergence of sums in the operator norm depends on a free analogue of the Burkholder- 
Gundy martingale inequalities. 
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A very useful tool in the study of stochastic integrals is the Ito formula. A free probability 
analogue of the Ito formula was developed in [BSa]. In terms of formal rules, it can be written 
as follows: 

atdt ■ btdt = atdt ■ btdWtCt = a t dWtb t ■ ctdt = 0, 
a t dW t b t -c t dW t d t = E(b t c t )a t d t dt. (9) 

Note that the rule in the second line is significantly different from the classical case. In terms 
of free stochastic integrals, the second rule can be written as follows: 

f 1 atdW t bt ■ f 1 c t dW t d t = f 1 ( f a T dW T b T )c t dW t d t 
Jo Jo Jo Jo 

+ f OtdWM fcrdWrdr) 

Jo Jo 



+ 



/ E(b t ct)a t dtdt. 
Jo 



(Compare Theorem 4.1.2. in MB Sal or Proposition 8 and Corollary 10 in (Aj.) Here is an 
illustration (a particular case of Proposition 4.3.2 in [BSa]): 
Let Wt be the free Brownian motion and define 

d(W?) := W?~ l dW t + W?- 2 {dW t )W t + ... + {dWtjW?- 1 . 

Then, 

(w a r = [ a d(wn 



[ a d{w?)+ f a V w?~ k - 1 - 2 (dw t ) w t k (dw t ) wl 

Jo Jo „^,~7i- o 



0<k+l<n-2 



and it follows that 

[»/2j-l 



r d(Wl l ) = (W a ) n - V in -2k- l)C k P W t n - 2k ~ 2 t k dt, 
Jo r - ^ 



(10) 



fc=0 

where C k are the Catalan numbers, 



An analogous formula for the classical Ito integral with respect to the Brownian motion B t 
is quite different: 

jf nBr\dB t ) = B n a - ^Z21 £ B r 2 dt. 
Below, we will use the free Ito formula in order to compute the moments of variable X t . 
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3. Free Stochastic Differential Equations 
3.1. Existence and uniqueness . A free stochastic differential equation (free SDE) 

dX t = a(X t )dt + b(X t )(dW t )c(X t ) (11) 
is a convenient shortcut notation for the following integral equation: 

X t = X + [ a{X T )dr+ f b(X T )dW T c(X T ). (12) 
Jo Jo 

We consider only equations with the coefficients that do not depend explicitly on time, and we 
will always assume that a(X t ), b(X t ), and c(X t ) are locally operator Lipschitz functions. (A 
function / : R — > C is called locally operator-Lipschitz, if it is a locally bounded, measurable 
function, and if for all A > 0, there is a constant Ka > 0, such that 

\\f(X)-f(Y)\\<K A \\X-Y\\ 

for all self-adjoint operators X and Y with the norm less than A. For example, all polynomials 
are locally operator-Lipschitz.) 

Equations (fTTT) and (PT2l) are particular cases of the following more general equations: 



m 

£< 

i=l 



dX t = a(X t )dt + > bi(X t )(dW t )ci(X t ) (13) 

and 



rt m rt 

X t = X + / a(X T )dT + Y] / bi(X T )dW T Ci(X T ). (14) 
Jo i=1 Jo 

Our results for the Cauchy transform of X t can be extended to this more general setting at 
the expense of more cumbersome notation. 

The existence of the solution of equation (fT4l) may fail for large t if the norm of the solution 
approaches infinity in finite time. However, for sufficiently small t > we have the following 
local existence result. (See Theorem 3.1 in MBSbl for a sufficient condition of the global 
existence in a simpler class of free SDEs, and Theorem 5.2.1 in (O) for an existence and 
uniqueness result in the case of classical SDEs). 

Theorem 3.1. Suppose that a\, bi, and Ci are locally operator Lipshitz functions and X is 
bounded in operator norm. Then, there exist to > and a family of operators Xt defined for 
all t G [0, to) an d bounded in operator norm, such that Xq = X, and Xt is a unique solution 
oftBfort<t . 

Proof: The proof proceeds by Picard's method of successive approximations. We will give 
the prove for the case m = 1. The general case is similar. Define X t = X, and 

X[ N+1) =X+ fa{xW)dr+ f \{xW)dW T c{xW). (15) 
Jo Jo 

We aim to show that this process converges for all sufficiently small t. For this, it is enough to 
show that for a sufficiently small to > and all t < to and N > 1, the following two claims 
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hold: (i) 

for some constant C and R, and (ii) 



x\ N) - xl N ~ i] 



A' 



(AO 



< c 



< A 



R 

ah 



for a constant A. 

Claim (ii) follows from claim (i) because (i) implies that 



X 



(AT) 



X 



<Cf(t) 



where / (i) := ^f* fc/2 is defiled for all t<l/R, monotonically increasing, differen- 

tiable at 0, and vanishes at zero. This implies that for all A > \\X\\ , there exists such to > 



that 



X 



(A) 



< A for all t < to- Moreover, this choice of to is independent of N. 



In order to prove (i), we proceed by induction. The case iV = 1 is special and can be easily 
verified separately. Assume that (i) and (ii) hold for x[ N ^ and x[ N ^ and let us prove that 
(i) holds for x[ N+l) . We write 



X 



(N+l) 



(AO 



< 



* r- 



(X^)-a(X^) 
b(xW)dW T c(xW) 



dr 



(16) 



The second term in this expression can be estimated by the following sum: 



2 \V2 

dr 



b(XW) - b (X^- 1 ))] dW T c(X^) 

+ / * b(X( N -V)dW T [c(X^) - c(Xi N -V) 
< 2V2 (J* \\b(xW) - b (x^- 1 )) || 2 \\c(XW) 

+2V2 (J* \\c(XW) - c (x^- 1 )) || 2 \\b(X^) 

where we used the free Burkholder-Gundy inequality (see Theorem 3.2.1 in HBSall ). 

By using the assumption that b and c are operator Lipschitz and claim (i), we see that this 
expression is bounded by 



2 \ V2 

dr 



pA 



T N dT 



1/2 



32K\A 2 C 



R 



N 



(N + l)\ 



N+l 



1/2 



A similar estimate can be obtained for the first part of (fl6l ). and by worsening a constant, we 
can obtain the following inequality: 

2 „ „ R N 



X (N+1) _ X (N) 



< UK\A 2 C 



(iV + 1)! 

{N+l) 



N+l 



provided that t < t . This shows that claim (i) holds for X^"' r±> provided that R > MK\A 2 
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Hence, the sequence X { t N) = [x[ N) - x[ N X) ) +. . . + (x[ l) - X^j +X is convergent in 
operator norm for every t < to- Let the limit be denoted by Xt- By using the free Burkholder- 
Gundy inequality, we can take limits on both sides of (fT5T > and check that X t is a solution of 

Next, suppose that X t and X[ are two different solutions of ([T2l for t < to- Let v (t) = 
\\Xt — X' t \\ . By using the assumption that the coefficients are operator Lipschitz and by using 
the free Burkholder-Gundy inequality, we obtain: 



v (t) < c\ J v (t) dr + c 2 v (t) 2 dr^j 

_ (Ft \ 1/2 

< {ciVk+C 2 ) U v(T) 2 dT 



1/2 



where c\ and C2 are certain positive constants that depend on Lipschitz constants. (The second 
inequality follows by the Cauchy-Schwarz inequality) By the Grownwall inequality (see [O], 
exercise 5.17 on p.80) it follows that v (t) 2 = for all t < to. Hence X t = X[ and we 
established the uniqueness of the solution. QED. 

3.2. Equations for the Cauchy Transform. 

Theorem 3.2. Assume that a, b, and c are locally operator-Lip schitz functions and let Xt be 
a solution of equation Alii bounded in operator norm for all t € [0, to)- Gt and gt denote 
the resolvent of Xt and the expectation of the resolvent, respectively, and let at = a(Xt), 
bt = b(Xt), and Ct = c(Xt). Then, for all t £ [0, to). 

^ = -E{a t G 2 t ) + E(btCtGt)E{btc t G 2 ). (17) 

Let us mention an important particular case, when the product be does not depend on X t . 
In this case, the equation simplifies to the following: 

^ = -E(atG 2 t) + [E{bc)] 2 gt d -^. (18) 

If we assume in addition that b = c = 1 and a(x) is a polynomial then (fT8T ) implies the free 
Fokker-Planck equation (O of Biane and Speicher. We will demonstrate this in Corollary 13.61 
below. 

In the proof of Theorem l3.2l we need the following lemma. 

Lemma 3.3. Let operators Hi and H 2 belong to the subalgebra W a which is generated by 
{W T } where r < a. Then 



E 



^ b T {dW T )c^j H x (J b T {dW T )c^j H 2 
E(c T H 1 b T )E(c T H 2 b T )dT. 
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This result follows if we write the integral as the limit of sums and use formula ([8]>. 
Proof of Theorem: For conciseness of the following formulas, let us use the following 
notation: 

/•t+At 

A = / a T dr, 

and 

/t+At 
b T (dW T )c T . 

Note that \\A\\ 2 = O(At) and ||5|| 2 = 0(y/At) for small At. 
By using the resolvent identity twice, we can write: 

Gt+At — Gt = -G t +At(A + B)G t 

= -G t AG t - G t BG t + G t+At (A + B)G t (A + B)G t 

Note that 

\\G t+At AG t AG t + G t+At AG t BG t + G t+At BG t AG t \\ 2 = o(At). 

In addition, 

||G , t +At -G t || 2 = 0(v / Al), 

which implies 

\\G t+At BG t BG t - G t BG t BG t \\ 2 = o(At) 

Hence, we can write 

E{G t+At - G t ) = E(-G t AG t - G t BG t + G t BG t BG t ) + o(At) (19) 
Next, we use the facts that 

-t+At 



that 

and that 

E 



pt+i\t 

j a T dr = atAt + o(At), 

/t+At 
G t b T {dW T )c T G t \ = 0, 



/ rt+ A t \ / rt+ A t \ 

gAJ b T (dW T )c T \G t U b T (dW T )c T \G t 



= (At)E(c t G t b t )E(c t G 2 t b t ) + o(At), 

where the latter holds because of Lemma [331 and the assumption that b t and q are Lipschitz. 
Hence, after taking the expectation in ( fT9l we obtain 

9t+At ~9t = At{-E(a t G 2 t ) + E(c t G t b t )E(c t G 2 b t )} + o(At), 

which is equivalent to the statement of the theorem. QED. 

In order to proceed further and obtain a differential equation on g t , we need to impose 
additional conditions on at, b t , and q which would allow us to eliminate expectations from 
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Proposition 3.4. Let Xt be the solution of equation riiil) , and G(t, z) and g(t, z) be its re- 
solvent and the expectation of the resolvent, respectively. Suppose that functions a and be are 
polynomials in one variable and that their degrees are not greater than k > 0. Then, 

dg d{ag) ^ (k - 1 - j)\ d^ +2 a(z) ■ 

Tt = k\ Jd^ E{X) m 

3=0 y 1 

3=0 '' K J 

3=0 y ' 

This equation is more useful than it might seem at the first sight. First of all, it is often 
possible to compute the expectations E(Xf)by using the Ito formula. Second, if these expec- 
tations are known, then the equation is a quasilinear PDE and the method of characteristics is 
applicable. 

Proof: Let f(x) be a polynomial. If we expand f(X)(X — z)^ 1 and f(X)(X — z)~ 2 in 
partial fractions and then take the expectations, we obtain the formulas: 

/(X), _ f(z) . , k ^(k-l-j)ldi +1 f(z) 



E{ x—l ] ~ E{ x—l ] + 2-* — *i — W 5 ^ ( } ' 



and 



3=0 



{X-zf 1 ' \X-zf' y X-z 

i-i)!^ +2 /( 

kl {dz)'J+ 2 



^ (fc -l-j)!ffl-+2 /(z) 

+ Z. ti E ( X3 )- 



3=0 

By using a(z) or 6(z)c(z) as /(z) it is easy to see that the statement of the proposition follows 
from Theorem |3.2j QED. 

Corollary 3.5. Suppose that a is a polynomial and that be = 1. Then, 
dg _ d(ag) [ dg 
dt dz dz 

3=0 V ' 



Corollary 3.6. Suppose that a is a polynomial with real coefficients, that b = c = 1, and that 
Xq is self-adjoint. Assume that the spectral distribution of Xt is absolutely continuous and 
bounded with the density p (x, t) . Then, at all points where dp/dx is defined, it is true that 

dp d , TT , 

-dt = -d-x {aP + P ' Hp} > 
where Hp is the Hilbert transform of p. 
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(This is the free Fokker-Planck equation (O of Biane and Speicher.) 

Proof of Corollary I3.6t Note that X t are self-adjoint for all t. Let us take the imaginary 
part on both sides of the formula in Corollary 13.51 and then pass to the limit y — > 0, where 
y := Imz. Assume that g(z) is analytic at z = x and therefore taking the limit commutes 
with operations of differentiation with respect to t and z. 

Since Xt is self-adjoint, therefore lm(E(Xj)) = 0. Hence, the formula in Corollary 13.51 
simplifies as follows: 

n tt=^l { -Yz lm[a9]+lm[9 ^z ]h (25) 
where we used the Stieltjes inversion formula. Note that 

lim Im(q— — ) = lim-flmgRe— — + Realm— — } 

y^o yij dz 1 y^o 1 dz y dz 1 

= n(p^(-Hp) - (Hp) A,) 



Similarly, 



lim -^-Jm[ag] = ir-^- [ap] 
J/->0 oz ox 



because Ima(x) = and Kea(x) = a(x). Hence, equation (l25l) simplifies to 

- = --[ap + p.Jp]. 

QED. 

3.3. Examples. In this section, we calculate explicit solutions in several particular cases. 
3.3.1. Ornstein- Uhlenbeck. 

Proposition 3.7. Suppose that Xt satisfies the equation of the free Ornstein-Uhlenbeck pro- 
cess: 

dX t = OXtdt + adWt. 
Suppose that Xq = 0. Then the spectral distribution of X t is the semicircle distribution 
supported on the interval Ig , where 
1) 



if 9 > 0, 
2) 

if0 = 0, and 
3) 



if 6 < 0. 



2rr 2 2rr 2 
^-(e^-l),+J^-(e^-l)], 



-2aVt, 2aVt] 



h 
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Hence, if 9 > 0, then the support of the distribution grows exponentially; if 9 = 0, then 
the support grows linearly, and if 9 < 0, the spectral distribution converges to the semicircle 
distribution supported on the interval [— aW2/\9\, <jy/2/\9\]. 

Proof: In this case at = 9X t , bt = \fo. Note that 

at = 6X t = 9{z + Gt 1 ). 

Hence, 

a t G 2 t = 9{zG 2 t + G t ) 

and 

E{atG 2 t ) = 9(zt^ + gt). 
Therefore, the differential equation for g t is 

£ + = (26) 

The initial condition Xq = corresponds to g(0,z) = —z , and we can solve this partial 
differential equation by using the method of characteristics (see pp. 9-19 in Q). 
Indeed the equations of characteristic curves are 

dt 

% = 0z-a 2 g, (28) 
d£ 

% = -0g. (29) 
By using d27l ), we can set £ = t. Then d29l) implies that 

and then we can solve d28l as 

z<t) = OP + ^ e -«. 
It follows that the initial point of a characteristic curve is given by equations: 

5(0) = A, 

,(0) = c + ^. 

On the other hand we can parameterize the initial condition of the PDE as follows: 

z(s) = s, g(s) = -l/s. 
Hence, we obtain the following paramterization for A and C : 
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Therefore, the equations of the characteristic surface are 



g(s,t) = --e- et , (30) 
s 

29 s 20 s 



From (1301 we have 

1 



ge et ' 



After we substitute this in (1311) and re-arrange, we obtain the following functional equation 

for g(t,z): 



g2 + a 2 (e 2et -1) 9+ a 2 (e 2et - 1) °' 
provided that 9^0. We can easily solve this quadratic equation for g. Note that by the Stielt- 
jes inversion formula the density of the corresponding distribution is given by the imaginary 
part of the Cauchy transform g. We can check that in our case this density corresponds to the 
density of the semicircle distribution. The radius of the semicircle distribution is 



29z 26 



2(T 2 

^r(e 2et ~ 1) 



if 9 > 0, and 




-2\0\t) 



if 9 < 0. This implies the statement of the proposition for 9 / 0. The case 9 = can be 
analyzed similarly. QED. 

3.3.2. Geometric Brownian Motion. Now let us consider the case when the coefficient b t 

1/2 

explicitly depends on X t . Namely, let at = 9X t , and b t = X t . 
In this example we deal with the equation 

dX t = 9X t dt + xl /2 (dW t )xl /2 , 

which is an analog of the classical equation for the "geometric Brownian motion", dxt = 

Oxtdt + xtdB t . 

Let us assume that Xq = I and use the free Ito formula to study the moments of the 
solution. Clearly, E(X t ) = e et . In order to calculate the second moment, we write 

d(Xf) = (X t + dX t ) 2 -X 2 

= (20X 2 + e 8t X t )dt + xf /2 (dW t )xl /2 + xl /2 (dW t )xf /2 , 

where we used the free Ito formula to calculate 

dX t dX t = x] l2 (dWt)X t (dWt)X] 12 

= E{X t )X t dt 

Let ht denote E{X 2 ). Then we have the following equation: 

^ = 29h t + 
dt 
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with the initial condition ho = 1. The solution is 

h t = (t + l)e 2et . 

Hence, the variance of the spectral distribution of X t is te 2et . The ratio of the standard 
deviation to the expectation of X t is 

In order to recover the entire spectral distribution, we use Theorem 13.21 and obtain the 
following result. 

Proposition 3.8. Suppose that Xt satisfies the following equation: 

dx t = ex t dt + xl /2 (dw t )xl /2 , 

and that Xq = I. Then, the expectation of the resolvent satisfies the following functional 
equation: 

z + g- 1 =e («-*ff)* j (32) 
where a = 9 — 1. The density of the spectral distribution of Xt is supported on the interval 

n(t) _( a _ n ( t )) t r 2 {t) (g- r2 (t))t] 



[ i + n {t) 'l + r 2 (t) 



where 



-i± v / r+4A 

n,2\t) = — 



2 

We can see from this proposition that the solution of the free SDE exists remains positive 
definite for all t > 0. 

If t — > oo, then r\^,{t) are asymptotically 1/t and —1 — \jt. Hence, as t — > oo the support 
of the solution becomes asymptotically close to 

l et ' J 

In particular, if 9 < 0, then both the lower and the upper bound of the spectral distribution 
shrink to zero exponentially fast, although at different rates (6 — 1 and 9). If 9 = 0, then the 
lower bound shrinks to zero exponentially and the upper bound grows linearly. If 9 € (0, 1), 
then the lower bound shrinks exponentially and the upper bound grows exponentially. If 
9 = 1, then the lower bound declines as (et)^ 1 and the upper bound grows exponentially. If 
9 > 1, then both the upper and lower bounds grow exponentially. 

Proof of Proposition 13.81 We have 

E{b* t G t b t ) = E{G t X t ) = 1 + zgt, 

and 

E(b*tG 2 bt) = E(G 2 X t ) = g t + z^. 

Hence the differential equation is 

d 9 at i d 9\ , n , \r , d 9t^ 
9t=-^ + z d - z ) + (l + z 9t )(9 t + z—), 

or 

£| + z{{e - 1) - zg)^ = -g((9 - 1) - zg). (33) 
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By assumption, the initial condition is g(0, z) = (1 — z) . 
The equations of characteristic curves are 

dt 



« 1. <34) 

dz 

= z ( d -l-zg), (35) 

^ = -g(0-l-zg). (36) 



From (1341 ) we can set £ = t. Then, if we divide (1361 ) by (1351 ). we obtain the following equation: 

dg_ = _9 

dz z 1 

which implies the following family of equations for the characteristic curves. 

A 

9 = —■ 

z 

If we substitute this in equation 051 ) and solve the resulting ODE, we find: 

z(t) = C e («-i-*)t. 

Hence, 

ff (t) = (A/CTJe-^- 1 -^. 
In particular, z(0) = C, g(0) = A/C. 

On the other hand, we can parameterize the initial condition of (l33l as follows: 

z{s) = s, g(s) = . 

1 — s 

This implies the following parameterization for A and C: 

s 



C = s, A 



1-s 

Hence, the characteristic surface is 



z(s,t) = S exp{(0-l--^-)i}, (37) 

1 — s 

g(s,t) = J—exp{-(0-l--^)t}. (38) 
1 — s 1 — s 



We can eliminate s from these equations: 



l + zg 

After we substitute this expression for s in (l37l) and re-arrange the terms, then we obtain the 
following equation: 

z + g- 1 =exp{(9-l-zg)t}. 

Let us denote 8 — 1 as a for simplicity of notation. Then, the functional equation for the 
Cauchy transform g(t, z) is as follows: 

z + g- 1 = e ( Q -^) i . (39) 
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If we take the differential of this equation, then we find that 

dz{\ + gte {a - z9)t ) = dg{g- 2 - zte^ 9 ^). 

The branch points of the function g(z) can be found from the equation dz/dg = 0. Hence, at 
the branch points, 

(a-zg)t _ 1 



g 2 zt' 

Substituting this into equation (|39l ), we obtain the following equation for the branch points: 

t(zg) 2 + t(zg) -1 = 0. 

Hence, 



-l ± JT+Ilt 
zg = \ '- = n, 2 (t). (40) 

Then (l39l > and (l40l imply that at the branch points, 

9i,2 = (l+ri, 2 (t)) e -( Q -^W)' 

and 

Zl „ - r J£$l p(«-n,2(t))« 
1,2 l + r li2 (t) e 

Finally, note that branch points of the Cauchy transform are bounds for the support of the 
spectral probability distribution. QED. 

3.3.3. Geometric Brownian Motion II. In our next example, we consider a different analog 
of the classical geometric Brownian motion equation, namely, the following free SDE: 

dX t = eX t dt + X t dWt + (dW t )X t . 

As in the previous example, assume that Xq = I, and note that E(X t ) = e et , the same as in 
the previous example. 

It is possible to write a PDE for the expectation of the resolvent in this example similar to 
equations in Theorem l3.2l and Proposition 13.41 However, it seems that it is difficult to find an 
explicit solution of this equation and recover the spectral distribution function of X t . 

Still, it is possible to see that the behavior of the solution is quite different from the behavior 
of the solution in the previous example by studying the variance of the solution. By using the 
free Ito formula, we can write: 

d(X?) = [29X? + X? + 2E(X t )X t + E(X?)]dt 
+X 2 dW t + 2X t {dW t )X t + {dW t )Xl 

Let ht denote E{X 2 ). Then we have the following ODE for h t : 

,U " 2(9 + l)ht + 2e 29t . 



dt 

The initial condition is ho = 1 and the solution is 

ht = 2e 2 ^ - e 2et 
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Hence the variance of X t is 2e (e — 1), and the ratio of the standard deviation to the ex- 
pectation is y / 2(e 2 ' — 1). This ratio grows exponentially fast with t, quite unlike the previous 
example, where this ratio equals \fi. 

3.3.4. Explosive equation. In our final example, we will consider an equation whose solution 
explodes in finite time. By this we mean that the norm of the solution becomes infinite in 
finite time. 

Proposition 3.9. Suppose that Xt satisfies the following equation: 

dX t = kX t (dW t )X t , 

and let the initial condition be Xq = al. Then the spectral distribution of Xt is defined for all 
t < (ak)~ 2 and it is supported on the interval: 

.(1-akVt) 2 (1 + aky/t) 2 , 



{l-a 2 k 2 t) 2 ' (l-a 2 k 2 t) 21 ' 

-i x 



I 

For t E (0, 1), the density of the spectral distribution of the operator a~ l X^ ak y T is given by 
the formula: 

m , V-(l-r)^ + 2(l+r)$-l 
f{0 = 2^~r • 

Proof of Proposition [32D We can compute 

E(G t b 2 ) = k[E(Gt 1 ) + 2z + z 2 g t ] 
= k[a + z + z 2 g t ], 

where we used the fact that 

E{G^) = E{X t ) -z = E(X ) -z = a-z. 

In addition, 

E(G 2 b 2 ) = k[l + 2z gt + z 2 ^]. 
Hence, the differential equation for g t is 

|| - k 2 (a + z + z 2 9)z 2 ^ = k 2 {a + z + z 2 g)(l + 2zg), 

and the initial condition is go = (a — z) . 

The equations for the characteristic curves are 

dt 

= 1. (41) 

d£ 

J = -k 2 (a + z + z 2 g)z 2 , (42) 
«? 

% = k 2 (a + z + z 2 g)(l + 2zg). (43) 
at, 
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From (|4TT) . we can set £ = t, and the equations for the characteristic curves in (z, g)-plane 
become: 

§ = -k\a + z + z 2 g)z\ (44) 

f| = fe2(a + z + ^)(l + 2^ 5 ). (45) 
as 



After dividing (1451 ) by (|44l . we obtain: 

(ig 1 + 2gz 
dz z 2 
The general solution of this equation is 

g (z) = -z- 1 + Cz~ 2 . (46) 

If we substitute this expression in (1441) . we obtain: 

Hence, 

*® = ^(a + ^-M " (47) 
By substituting this in (l46l >. we obtain: 

5 (t) = -(a + C)£; 2 r - A + C((o + C)& 2 t + A) 2 . (48) 

In particular, if i = 0, then 

z(0) = (49) 
5 (0) = -A + C.4 2 . 

On the other hand, the initial condition is <?(z) = (a — which we can parameterize 
as follows: 

z{s) = s, g(s) = (a-s)- 1 . (50) 
Comparing (|49l and (l50l . we obtain the following parameterization for A and C: 

,1 i.r "* 



s a — s 

We substitute these expressions in ( |4"7T ) and d48l ) and obtain: 

*(*>*) = -T— n (51) 

</(*,*) = r ^fc 4 t 2 + ^±^fc 2 t + — . (52) 
(a — s)- 3 [a — s) z a — s 

We are going to eliminate s from the pair of equations (TSTI ) and (l52l . For this reason, we write 
(l52l as follows: 

s(t, .) = -^(^fc 2 t + ±)(-2-A*t + 1), 
a — s a — s s a — s 

and then we substitute ( f5TT > and obtain: 

5 (t,s) = -(^-fc 2 t + 1), 
a — s z a — s 
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or 

,3 



S 

By using (|5~T1) again, we note that 



Z5 = /c 2 i + l. (53) 



a 3 , o „ a a 
k 2 t + 1 = + 1. 



a — s z s 

Hence, (l53l can be re- written as follows: 

a 

(a — s)zg = ( — h l)s — a, 

z 

and, therefore, 



s = a: 



a ' 



and 



1 + zg + § 

After substituting these expressions in equation d5TT >. we obtain: 



I = ( z + z 2 5 + a)k 2 t + - + ' 



z a z + z 2 g 

After re-arranging the terms and dividing by z, we get the following equation: 

z 1 

k 2 tz 3 g 2 + ( 1 + (a + 2z)zk 2 t)g -\ V (z + a)k 2 t = 0. 

a a 

This functional equation for g(z, t) is quadratic and therefore it is easily solvable. 

In particular, the branch points of gt{z) are the zeros of the discriminant of this equation, 
which can be computed as 

D = (k 2 at - -) 2 z 2 - 2(k 2 at + -)z + 1. 
a a 

Therefore, the branch points are 

(1 ± akVt) 2 
Z± ~ a {l-a 2 k 2 tf 

Note that as t approaches (ak)~ 2 , the branch points approach a/4 and oo. 

It follows that for t < (ak)~ 2 , the spectral distribution of X t is supported on the interval 
z+] and in this region it has the density 



/(x)da; = -- * a ^ da:. 

a 2 / nk z a z t{x/a) 6 

If we use variables r = k 2 a 2 t, and £ = x/a, then we can write this density as 

^-(1 " t) 2 { 2 + 2(1 + r){ - 1 



QED. 
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